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Abstract. Unitary transformations and density matrices are central objects in 
quantum physics and various tasks require to introduce them in a parameterized 
form. In the present article we present a parameterization of the unitary group 
U(d) of arbitrary dimension d which is constructed in a composite way. We show 
explicitly how any element oilA(d) can be composed of matrix exponential functions of 
generalized anti-symmetric er-matrices and one-dimensional projectors. The specific 
form makes it considerably easy to identify and discard redundant parameters in 
several cases. In this way, redundancy-free density matrices of arbitrary rank k can be 
formulated. Our construction can also be used to derive an orthonormal basis of any 
fc-dimensional subspaces of C d with the minimal number of parameters. As an example 
it is shown that this feature leads to a significant reduction of parameters in the case 
of investigating distillability of quantum states via lower bounds of an entanglement 
measure (the m-concurrence) . 



A composite parameterization of unitary groups, density matrices and subspaces 2 
1. Introduction 

In quantum information many quantities or properties of quantum systems are 
related to optimization problems that necessitate to vary over the set of all unitary 
transformations, density matrices or subspaces. Such problems arise, for instance, 
in the detection and quantification of entanglement (see [H EJ [3] and references 
therein), the properties of quantum states with respect to Bell inequalities [U E] or 
the question of distillability [7]. In general, such a variation can only be done 
efficiently by parameterizing the object of interest. In the present work the focus is 
on the parameterization of the unitary group U(d). Density matrices and orthonormal 
subspaces can then easily be formulated in terms of this representation. It should 
be emphasized that parameterizations are in general equivalent to each other and one 
should not be misled to conclude that one parameterization is better than the other. 
However, depending on the problem a certain parameterization can be advantageous by 
providing more insightful results or by leading to a reduction of the number of involved 
parameters. The simplest parameterization of U(d) is the canonical parameterization 
which is given by U = exp(iFT), wherein the hermitian matrix H can be composed 
of a real-valued linear combination of the d 2 — 1 generalized Gell-Mann matrices and 
the unity 1^. Despite its simple form, this parameterization has the disadvantage that 
there is no way how redundant parameters, which appear for some tasks can be removed 
beforehand (except for one parameter for the special unitary group SU(d)). Thus the 
number of parameters is always d 2 or d 2 — 1, respectively, independent from the given 
problem. A different parameterization of U(d) and SU(d) in terms of generalized Euler 
angle was introduced by Tilma and Sudershan [HI E]- It has the advantage that it 
allows to eliminate redundant 'unphysical' global phases in several cases. A further 
parameterization of U(d) was recently found by Jarlskog [TUl [TTj . Here, the unitary 
matrices are formulated in a recursive way, meaning that the elements of U(d) are 
expressed in terms of the elements of Hid — 1) and a further unitary matrix containing 
the parameters which enables the extension to d. This parameterization also enables 
one to remove invariant phases, however, this is a nontrivial task as shown in [T2] . 

In the present work we introduce a parameterization (see Sec. 2) which is ideal for 
the formulation of density matrices (see Sec. 3) and orthonormal subspaces (see Sec. 4) 
since all redundancies can be easily identified and removed beforehand without having 
to consider fixed dimensions or special cases explicitly. Due to its concise notation and 
formulation in terms of matrix exponential functions of dyadic vector products, the 
parameterization can be easily implemented in computational programs. Moreover, the 
parameters have simple ranges; when gathered in a d x d matrix in a particular order, 
they permit an insightful interpretation of numerical computations. Concrete numerical 
examples, in which our parameterization reduces the number of involved variables are 
given in Sec. 5. Here, we state which parameters can be discarded with respect to 
optimizing lower bounds of an entanglement measure (the m- concurrence), as well as 
for clarifying distillability of quantum states. 
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2. Composite parameterization of the unitary group U(d) 

Consider unitary operators U acting on a Hilbert space H = C d with d > 2 spanned by 
the orthonormal basis {|1} , . . . , \d)}. For any U G U(d) there exist d 2 real values \ m , n 
with m,n E {h . . . , d} and \ m , n £ [0, 2tc] for m > n and \ m ,n £ [0, f ] for m < n such 
that U equals | 



c 



II EI exp ( iP ^n, m ) exp (icr mjn A mi „) 



_m=l \n=m+l 

Here, the are one-dimensional projectors 

Pi = 10 CI 

and o- m , n are the generalized anti-symmetric cr-matrices 

o- m , n = —i I'm) (n\ + i \n) (m\ 



Y[exp(iPi\i t i) 



(1) 



(2) 



(3) 



with 1 < m < n < d. 

Before we prove that any unitary operator of U(d) can be written in the form (pQ) let 
us comment on the concept behind it: Unitary transformations have the characteristic 
trait that they map orthonormal bases onto orthonormal bases, i.e. for a given set of 
vectors , . . . , \d)} forming a orthonormal basis of a Hilbert space C d the unitarily 
transformed set {|1') , . . . , \d')} = {U |1) , . . . , U \d)} is also orthonormal. We are 
interested in the most general unitary operations transforming {|1} , . . . , \d)} into any 
arbitrary orthonormal basis , . . . , \d')} , hence our object is the unitary group of 
dimension d. We know that for dimension d the unitary group U(d) is a <i 2 -parameter 
group. This means that in order to cover all unitary transformations we must find a 
transformation that contains at least d 2 parameters (if a parameterization involves more 
than d 2 parameters it also contains redundancies which is undesirable). 
Regarding the basis {|1) , . . . , \d)}, we first find that d parameters can be embedded 
in global phases for each vector, i.e. {e iai |1) , . . . , e lotd \d)}. This corresponds to the 
product FJ^ =1 exp(iP/A; i ;) of our parameterization ([T|). The notion behind the left part 
of ([!]) is to pairwise embed two parameters in the distinctive transformations 



^-m,n exp {iPn^n,m) GXp iio~ m ^ n A m>n ) 



(4) 



containing two parameters \ n , m and \ m ,n- These transformations perform the following 
operations: The term exp {ia mtn \ myn ) generates a rotation in the two-dimensional 
subspace spanned by the vectors \m) and |n), while exp {iP n X n)m ) adds a relative phase 
between the vector components of the rotated vectors. Note that U (2) is actually a 4- 
parameter group, however, as will be proven later, the neglected two parameters would 
only lead to redundancies. In order to parameterize U(d) all terms A m n have to be 



| The sequence of the product is YiiLi A-i = A\ ■ A?, - ■ ■ A 



N 
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taken into account. There are 



d(d— l)/2 ways how two vectors of {|1) , . . . , \d)} 



can be combined, which corresponds to the d(d — l)/2 generalized anti-symmetric a- 
matrices. One opportunity to involve all terms A mj „ is given by building the product 
(l~In=m+i ^-m,nj , which is clearly not unique. We now have embedded further 
2 x d(d — l)/2 = d 2 — d parameters, which in sum with the d global phases gives d 2 
parameters \ m ,n m total. For convenience, these parameter A mi „ are gathered in a d x d 
'parameterization' matrix 



/ Ai,i 



(5) 



\ A 



A, 



dA 



In this convention the diagonal entries represent global phase transformations, the upper 
right are related to rotations, while the lower left are relative phases (with respect to 
the basis {|1) , . . . , \d)}). 

We have thus constructed a d 2 parameter set of unitary operators. It remains to prove 
that this construction covers the entire unitary group U(d), which is equivalent to 
showing that for any U there exists a Uc such that U C U = 1. 
Proof: 



E 



i=l a i,m a i,n 



Let U 
yd 



r,s=l fl »V 



\r) (s\ be an arbitrary unitary operator, i.e. Ylt=i a m,i a n,i 



5 mn . The conjugate transpose of Uc is 



n ex p(- ip d+i-i^d+i-i,d+i-h 



.1=1 



d-l 



II IK 

m=l \n=l 



—m,d+l—n 



(6) 



implying that A\ 2 acts first on U. For U' = A.\ 2 U = J2ts=i a 'r,s \ r ) ( s \ this leads to 



"l,s 



cos(Ai )2 )ai jS - e ^ sin(Ai i2 )a 2 ,5 , 



sin(Ai j2 )ai )S + e cos(Ai )2 )a 2jS . 



(7) 
(8) 



For d > 2, all other components are unchanged, i.e. a' r s = a r>s for r > 2. As can easily 
be confirmed, a' 21 can always be made zero via certain parameters A ij2 and A 2j i. If 
ai ; i and a 2ji both are zero, both parameters Ai )2 and A 2ji can be chosen freely; if only 
= we choose Ai j2 = | and for a 2i i = we take Ai j2 = 0. If none of both is zero 
then a' 2 x vanishes for A 2ji and Ai ;2 obeying 



arg(e lX2 ' 1 a 2 ,i) = arg(-a hl ) , 

l°2,l| 



tan(Ai j2 ) = 



°i,i 



(9) 
(10) 



which is achievable in any case with A 2j i G [0, 27r] and a Ai j2 G [0, |]. In the same 
way the component a'l : of U" = A} 3 C/' can be made zero for d > 2. By proceeding 



A composite parameterization of unitary groups, density matrices and subspaces 



5 



in this way for d > 2, we can attain a r>s = for all components with r > s via 
Ut=i (iXLi Ad-wm-n) ■ Then, if a rA = for all r > 1 it follows |a M | = 1 due to 



unitarity (J2 i= i a *,i a i,i = I) of U P = f] 



m=l 



FT" A f 

lln=l Jv d-m,d+l 



[/. Moreover, since 



C/p also satisfies J2t=i a i,i a i,i = 1 we can infer that a 1)T = for all r > 1. When taking 
into account the unitarity constraints Ylt=i a m,i a n,i " x " 
and columns we conclude that U p has the diagonal form 



Ei=l «*,m a ^™ = S mn for a11 rows 



d 

E 

r=l 



(11) 



with a rr obeying |a r>r | = 1, which is a complex number of magnitude 1, i.e. a r>r = e %OLr . 
The choice \ r = a r then yields U ] C U = Y\f =l exp(-iP d+1 _i\ d+1 _i 4+1 _i) U P = 1, 
which was to be proven. 



3. Parameterization of density matrices with rank k 

Now we use our parameterization of U(d) to formulate density matrices. Any density 
matrix p acting on H = <C d can be written in the form 

k 

p = Y^Pn\^n)(^n\ (12) 

71=1 

with p n > and Ylt=iPn = 1 where k < d is the rank of p and {l^i) , . . . , |^fc)} are 
orthonormal vectors. For k — 1, i.e. pure states we only have one p\ = 1. Without loss 
of generality, the coefficients {p n } of an arbitrary state of rank k > 1 or smaller can be 
expressed by means of k — 1 real parameters 9i G [0, |] as 

Pl = cos 2 6 l (13) 

n-l 

p n = cos 2 9 n Y[ sin 2 9i Vn : 1 < n < k (14) 
i=i 

fc-i 

^ = n s ^ 2 ^ • (1 5 ) 

Any desired set of k orthonormal vectors , . . . , |^fe)} can be constructed out of 

{|1) , . . . , \k)} by applying Uc- Thus any p with rank k can be parameterized via 

k 

P = 5>nt/c |n> (n| f/J • (16) 

n=l 

A parameterization of density matrices via diagonal elements and unitary transforma- 
tions can in principle be achieved with any parameterization of the unitary group U (d) 
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(see for instance [13], [HI [15]). However, the composite form of Uc makes it easy to iden- 
tify and eliminate all redundant parameters. The first observation is that the diagonal 
entries A„ jn are redundant, since they cancel out in the outer product, i.e. 



1=1 



\n) (n| 



J^exp(-zP d+ i_iA d . 



+l-l,d+l-lj 



1=1 



\n) (n| Vn G {1, .., d} . 

(17) 



The second observation is that for density matrices of rank k < d we can eliminate 
further parameters when the composite parameterization is used. As can easily be seen 
p is independent of all parameters \ m ,n where both m > k and n > k. Thus, it suffices 
to utilize 



P = ^PnUcD \n) (n\ Ul 



CD 



n=l 



with 



U C D = H II eXP ( iP n^n,m) exp {iCT m ^ m ,n) 



(19) 



m=l \n=m+l 



where the index m is only running from 1 to A; instead of 1 to d— 1. Using the introduced 
representation in terms of the parameterization matrix, this is Uc with 



/ Ai, 2 
A24 

Afc+1,1 ■ 

V x d,i ■ 



Ai 



fc+i 



Afc^+i 
Ajfc+i fe 



A 



d,k 







Afe,d 


J 



d — k 



(20) 



d-k 



Consequently, any density matrix of rank k or smaller acting on C d can be expressed 
with a maximum of 2dk — k 2 — 1 parameters (k — 1 parameters 0$ and fc(2d — k — 1) 
parameters A m „). The extremal scenarios are pure states with 2(d— 1) parameters and 
full rank density matrices with d 2 — 1 parameters. 



4. Parameterization of /c-dimensional subspaces 

In the previous section we have shown that any set of k orthonormal vectors can be 
constructed (up to global phases) with k(2d — k — 1) parameters. Now, we prove that 
even less, namely 2k(d — k) parameters are necessary to construct an orthonormal basis 
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of an arbitrary k- dimensional subspace of C d . Consider a general subspace spanned by k 
orthonormal vectors , ■ ■ ■ , \ ^k)} in a d- dimensional complex Hilbert space % = C d . 
Let a m , n be the coefficients of \^f n ) in the complete basis , . . . , \d)} of % = C d , i.e. 
\^n) = Em=i a m,n \ m ) ■ A different basis of the same subspace then is of course given 
by the orthonormal set of vectors 

|* fc - 2 ) , cos(A fe _ lifc ) |* fc _!) - e (^.*-i) sin(A fe _ lifc ) |tt fc ) , (21) 
sin(A fc _ 1)fc ) |* fc _!> + e ^ fc ' k -^ cos(A fc _i,fe) |* fe ) } , (22) 

which is {A' fc _i, fc , ■ ■ • , A' fc _ lifc |* fe )} where A' k _ 

-l k * s gi ven as i n Q but with Pi — 
|*/) and a m>n = -i \^ m ) (^„| + % |* n ) (* m |. The (fc - l)t/i vector of this set is 



cos(A fc -i,fe) - < 

d 



(iXk,k-i) 



sin(A fc _i )fc ) |^ fc ) 



^2 (cos(A fe _i jfc )a miA ._i -e (iAfc ^- l) sin(A fe _i ifc )a mifc ) \m) 



(23) 
(24) 



m=l 



whose k£/i coefficient cos(\k-i t k) a k.k-i ~ e^ Xk - k - 1 ' > sin(Afc_i./ c )afc i fc can be made zero 
analogously to (l7l)-(fTUl). In this way, with {n n =i Ki,k l^i) 5 • • • 5 lln=i ^4 jfc l^fc)} a basis 
of the subspace can be obtained where all kth coefficients of all vectors except the 
last one are zero. This can then be repeated for all rows < k, such that with 

{U k m =2 Un=i A 'n, m , ■ ■ ■ , llm=2 Un=i K, m I**)} we arrive at an orthonormal basis, 
let us say {\^[) , . . . ,\^' k )}, of the same subspace where all coefficients a m<n = for all 
m with n < m < k. Such a set of vectors however, can be constructed (up to global 
phases) from {|1) , . . . , \k)} via {Ucs |1) > • • • > Ucs \k)} where 



U C s ] [ ) [ exp {iP n X nt m) exp (ia m>n \ min ) 



(25) 



m=l \n=k+l 



Here, only 2k(d—k) parameters A m n are involved since the index n of the second product 
Y[ lies within k < n < d. In terms of the parameterization matrix, this is Uc with 



/ 



Afc+1,1 



Ai, 



fe+i 



Afc^+i ■ 
Afc+i,fe 

Awfc 



Afc,d 


/ 



d — k 



(26) 



d-k 



Proof: 

The proof is completely analogous to the one in section 2 and will therefore not be 
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given in detail. One only has to take an arbitrary set of vectors {\^[) , . . . , l^)} 
with the mentioned properties, for which via {Uq S |^} , . . . , Uq S |^}} = 
{J2m=i a m,i \m) , . . . , J2 m =i a m,k \m)} it is possible to attain \a m>n \ = 5 m , n for all 
m G [1, . . . , d] and n G [1, . . . , k] - QED. 



5. Application: Optimization of lower bounds of entanglement measures 

The m-concurrence introduced in Ref. [TB] constitutes a building block of entanglement 
measures for multipartite systems of arbitrary dimension (see also Ref. [17]) with 
simple computable lower bounds. These bounds are not invariant under local unitary 
transformations. Optimization procedures can be realized efficiently via the previously 
introduced parameterization of the unitary group. 



5.1. Optimal lower bounds 

First we introduce the m-concurrence for bipartite quantum systems in C d £§> C d . The 
linear entropy Sl(p) := ^rx(l — Tr(p 2 )) of the reduced density matrices Pa/b '■= T^b/a{p) 
of a bipartite pure state \ip) G C d ® C d can be expressed as 



^J^-Sl(Pa) = E ^)^\^ A u®^ B M)(MT^k A u®^k B ,i^ (27) 

k A <lA k B <l B 

=-<%M)Q>\)- (28) 



Via a convex roof extension one can generalize the m-concurrence C^J^ip) to mixed 
states 

c M--= /^^PiCimm) (29) 

{pi,m)} 

= .in( u & E E ^(l^*)<V'<kfcx,U®^fl,iB(lV't)(V'i|)Vfc A ,i A ®<Tfc Bj j fl ). 



i Ua<Ia k B <l 



As the infimum of this expression cannot straightforwardly be computed, it is of great 
importance to find strong lower bounds. Those can be obtained with 



C(P)> E E inf E^'^^l 

0~k A ,l A ® cr k B ,l B {\i ! i) (lpi\)*0-k A ,lA ® 0-k B ,l B ) 

,,,,,, xPiAWi/S 
k A <iA k B <l B I 

by exploiting that the individual infima are known (see Ref. [TSl [T5] ) 
with 

x k A U,k B ,i B ■= msx[2max[{x\ A)U)kB)lg }] - E X L,U,<W S > °1 ( 31 ) 
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where {x\ 1a kg 1b } are the square roots of the eigenvalues of 

P °~k A ,i A ® crfc s ,« s P* <?k A ,i A ® crfc s ,/ s ■ (32) 
In summary we have obtained the following bound 

CI(P) >EE X Li A ,k B , lB =■ B 2 (p) . (33) 

k A <l A k B <l B 

This bound is not invariant under local unitaries and thus can be optimized with 
an appropriate choice U A <8> Ub- Consequently, the optimal lower bound for the m- 
concurrence for bipartite systems is given by 



BlpM ■= max 



E E 



k A ,l A ,k B ,l B 

,k A <l A k B <l B 



(34) 



where the eigenvalues of p o~ kA ^ A ®a kB ^ B p* c kA ^ A ®a kB ^ B are replaced by the eigenvalues 
of 

U A ®U B pU ] A ® U B a kAU g> o- kB ,l B U* A <g> ^bP #c/ I ® U B°k A ,i A ® (35) 

and the maximum is taken over all £7a, G Ll(d). In many cases, this problem can be 
solved numerically by parameterzing Ua, Ub and utilizing numerical methods such as 
'Nelder-Mead' [20], 'simulated annealing' [21J or 'differential evolution' [22]. In order to 
remove redundant parameters we exploit that 



EV \U A <g> U B pU\ ® U j B a kAilA ® <Tjfc B ,, B t^ ® C/bP*C/1 ® C/Jcr fcA; ^ ® a feBiJB J (36) 
=EV (p (u\a kAM U* A ) (S> P* (^x^^) ® (37) 

where EV stands for the set of eigenvalues. If we now insert U\ and U B in parameterized 
form (JTJ), i.e. C7^ = ?7c,a an d ZTg = Uc,b, and let them act on the cr-matrices, it is easy 
to see that all diagonal entries \ m ,m of both unitaries cancel out. Consequently, we 
can reduce fl34l) to a 2(d 2 — d) dimensional global optimization problem. (If we would 
insert U A = Uc, A and Ub = Uq,b, only one diagonal entry \ m ^ m of each Uc could 
generally be set zero.) An illustrative example, where the bounds of the m-concurrence 
are optimized in this way using the Nelder-Mead method [20] for convex combinations 
of two mutually orthogonal maximally entangled states and uncolored noise is given in 
Fig. [TJ With optimization the bounds are greater than zero for all states detected by 
the partial transposition criterion (PPT). 

For multipartite systems we can use the same principle (for a detailed introduction 
see Refs. [TH E]). First we introduce the set of all bipartitions B = {(ck|/3)} of a 
given n-partite system. Here, a denotes a union of subsystems in the first part of 
the bipartition and (5 is its complement. The dimensions of the corresponding complex 
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(a) 



a 




Figure 1. Contour plots of the lower bounds B(p) and B opt {p) of the normalized 
m-concurrence c n^T^Wn ^ or t ^ ie set °^ states P ~ a l v f'i}(*i| + /9 1*2} + 
1 ~q~ /3 1 constructed with two mutually orthogonal maximally entangled states = 
■j= (|11) + \22) + |33» and |* 2 ) = ^ (|12) + |23) + |31» of C 3 ® C 3 (bipartite qutrit 
system) combined with uncolored noise 1 . All values of the parameters a and /3 within 
the green triangle correspond to density matrices (positive semidefmite) . Density 
matrices within the blue ellipse are positive under partial transposition (PPT criterion). 
In the grey shaded areas the bounds of the m-concurrence are nonzero. The grey scale 
corresponds to the value of the bounds lying in the range of and 1. The shades of grey 
are related to an increment of 0.2 starting from 0. The left picture (a) illustrates the 
contour plot of the bounds computed according to equation (|33[) without optimization. 
The right picture (b) illustrates the contour plot of the numerically optimized bounds 
using our parameterization for U{2>) and the Nelder-Mead method. Without 
optimization (a), the lower bounds are zero for some states with negative partial 
transposition. After the numerical optimization (b), the lower bounds of all states 
that are negative under partial transposition are greater than zero. As can also be 
seen, the shapes of the contours of (b) differ considerably from (a), and indicate an 
improvement of the lower bounds. 



Hilbert spaces 7i a and Hp are d a and dp, respectively. In this way, the general definition 
of the m-concurrence also valid for multipartite systems is 

= , E ft £ E E Tr (i^> w<u ® <i B ) 

i a ka<Ia kb<Ib 

where {c^f l } and {a^ i B } are the generalized anti-symmetric a- matrices defined with 
respect to the bipartition (a|/3), i.e. acting on % a and Hp, respectively. Consequently, 
according to (J2HD - flM} 



^(P)>E E E (<U,Wb) 2 > 08) 



B kA<lA ks<lB 
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where X^ 1a kg lg is defined as in f[3T]) but with the eigenvalues of 

aB _ aB * aB „ aB /nc\\ 

P °kl,u ® a ks,i B P °k A ,i A ® • (39) 

By proceeding analogously to - ( 157)) . this lower bound can be optimized by replacing 
p by U A ®U B pU A ®U B f° r eacri bipartition. In this sense, the optimal lower bound 
is given by 

(40) 

over all XJf e U{d a ) and L/]f G U(dp). 
5.2. Distillation 

If entanglement is used as a resource it is often required that the system is in a 
pure maximally entangled state. States that can be transformed into such maximally 
entangled states via LOCC (local operations and classical communication) are called 
distillable. It was proven by Horodecki et al. in [6] that all distillable states have an 
entangled C 2 <8> C 2 subsystem (see also Ref. [23] )• As the bounds for the m-concurrence 
are exact in those systems it suffices to optimize one of the terms in the sum of (1341) 
to investigate distillability since all terms are local-unitarily related. Thus, a state p is 
distillable only if 

maxJT 2 212 > , (41) 

where the maximum is taken over all U A £ U{d) and U B EU(d). According to (1571) . the 
value of Ai 2,i,2 is a function of the eigenvalues of 

P (u\a lt2 U* A ) ® {yla^U^j p* (Ula li2 U A ) ® {U T B a^ 2 U B ) . (42) 

These eigenvalues, however, are completely determined by the 2x2 dimensional subspace 
spanned by the tensorproducts of the vectors U A |1) and U\ |2) with U B |1) and U B |2). 
Consequently, for W A and U B we can take U' A = Ucs,a and U B = Ucs,b defined as in (1251) 
with k = 2, where each transformation only depends on 4d — 8 parameters X m ,n- In this 
way, the question of whether a state p is distillable or not, can efficiently be clarified via 
a (numerical) optimization algorithm with a reduced number of parameters. (Note that 
in contrast to a naive parameterization of U(d), the number of parameters is linear in d 
instead of quadratic. As distillability is generally studied for high-dimensional n-copy 
states, i.e. p® n , this reduction is of great importance for numerical tractability.) 

Multipartite systems can be treated equivalently for each bipartition. Hence, for a 
fixed bipartition 4(d a + dp) — 16 variables have to be optimized. However, in this case, 
if the unitary operators related to a fixed bipartition are not locally implementable with 
respect to the subsystems, also unlockable bound entanglement can appear (see also 
Refs. [2H ES])- The set of locally distillable states can be determined by restricting all 
transformations to unitaries which are local with respect to all subsystems. 
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6. Summary 

In this paper we introduced a parameterization of the unitary group U(d) which beside 
its simple form has the advantage that redundancies can easily be identified and re- 
moved in several cases. The efficiency for the representation of density matrices of rank 
k and k- dimensional subspaces in C d was shown. We found that only 2dk — k 2 — 1 real 
parameters are necessary to parameterize density matrices of rank k since we were able 
to discard all irrelevant parameters related to transformations in uninvolved subspaces 
and invariant phases beforehand. For the construction of an orthonormal basis of any 
fc-dimensional subspace of C d even less parameters are needed, namely 2k(d — k), due 
to the unitary equivalence of basis vectors within the subspace. 

Furthermore, examples of the usefulness of the parameterization with respect to a mul- 
tipartite entanglement measure (referred to as the m-concurrence) and distillability 
were given. We described how lower bounds can be derived for this entanglement mea- 
sure. The bounds were obtained by the observation that the linear entropy of reduced 
density matrices can be rewritten by an operator sum, where each operator acts in a 
two-dimensional subspace. We showed how these bounds can be optimized and how 
invariant parameters can be removed when our parameterization is utilized. Finally, we 
revealed that further parameters can be discarded if only the distillability of quantum 
states is of interest. We found that only a number of variables linear in the dimensions 
of the subsystems has to be optimized to solve this problem. 

We believe that the parameterization presented in this paper is advantageous with re- 
spect to the tractability of a variety of high-dimensional optimization procedures in 
quantum information as well as for other problems in quantum theory. 
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